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We derive the general formula giving the Berry's Phase for an arbitrary spin, having both 
magnetic-dipole and electric-quadrupole couplings with external time-dependent fields. We assume 
that the "effective" electric and magnetic fields remain orthogonal during the quantum cycles. This 
mild restriction has many advantages. It provides simple symmetries leading to useful selection 
rules. Furthermore, the Ifamiltonian parameter space coincides with the density matrix space for 
a spin 5 = 1. This implies a mathematical identity between Berry's phase and the Aharonov- 
Anandan phase, which is lost for S > 1. We have found, indeed, that new physical features of 
Berry's phases emerge for integer spins > 2. We provide explicit numerical results of Berry's 
phases for S = 2,3, 4. For any spin, one finds easily well-defined regions of the parameter space 
where instantaneous eigenstates remain widely separated. We have decided to stay within them 
and not do deal with non-Abelian Berry's phases. We present a thorough and precise analysis of 
the non-adiabatic corrections with separate treatment for periodic and non-periodic Hamiltonian 
parameters. In both cases, the accuracy for satisfying the adiabatic approximation within a definite 
time interval is considerably improved if one chooses for the time derivatives of the parameters a 
time-dependence having a Blackman pulse shape. This has the effect of taming the non-adiabatic 
oscillation corrections which could be generated by a linear ramping. For realistic experimental 
conditions, the non-adibatic corrections can be kept below the 0.1 % level. For a class of quan- 
tum cycles, involving as sole periodic parameter the precession angle of the electric field around 
the magnetic field, the corrections odd upon the reversal of the associated rotation velocity can be 
cancelled exactly if the quadrupole to dipole coupling is chosen appropriately ( "magic values" ) . The 
even ones are eliminated by taking the difference of the Berry Phases generated by two "mirror" 
cycles. We end by a possible application of the theoretical tools developed in the present paper. 
We propose a way to perform an holonomic entanglement of A*' non-correlated one-half spins by 
performing adiabatic cycles governed by a Hamiltonian given as a non-linear function of the total 
spin operator S, defined as the sum of the A*' individual spin operators. The basic idea behind this 
proposal is the mathematical fact that any non-correlated states can be expanded into eigenstates 
of and S^.The same eigenvalues appear several times in the decomposition when N > 2 but all 
these states differ by their symmetry properties under the A'-spin permutations. The case N = 4 
and Sz = 1 is treated explicitly and a maximum entanglement is achieved. 

PACS numbers: 03.65.Vf, 03.67.Bg, 03.75.Dg, 37.25.-|-k 
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I. INTRODUCTION 

About twenty five years ago, Geometry made a new 
entrance in Quantum Mechanics with the discovery of 
geometric phases These developments hinge upon 

the simple fact that all the physical information rela- 
tive to an isolated system described by a pure quantum 
state l^") is contained in the density matrix p — 
This mathematical object has two important properties: 

i) it is invariant upon the abelian gauge transformation 
1^") — ^ exp(ix) IvP) where x is an arbitrary real number, 

ii) it satisfies the non-linear relation = p. The den- 
sity matrix space E(yo) has clearly a non-trivial topology. 
A geometric phase is acquired when the quantum sys- 
tem performs a time evolution along a closed circuit on 
E{p), and satisfies at every instant t the so-called "parallel 
transport" condition = 0. The above terminol- 
ogy reflects the fact that the set of quantum vector states 
associated with a given density matrix p can be viewed as 
the "fibre" of a linear "fibre bundle" constructed above 
the "base space" E(p). 

Our work deals mainly with adiabatic quantum cycles 
performed within a definite time interval < t < T. 
They are generated by an Hamiltonian depending on a 
set of parameters X(t), assumed to be a system of co- 
ordinates for a differential manifold M. In this way, an 
adiabatic quantum cycle generates a mapping of a closed 
circuit drawn upon the parameter space M onto a closed 
loop upon the density matrix space E(p). The Berry 
phases can then be viewed as the geometric phases asso- 
ciated with this particular class of quantum cycles. 

At this point, a question arises naturally; how can 
one extract the geometric phase since p is phase in- 
dependent? The answer lies within the Superposition 
Principle of Quantum Mechanics: any linear combina- 
tion of two quantum states and [^'2) relative to 
a given quantum system, |^'i2) = ci|^'i) -I- C2|^'2) is 



an accessible state for the system. In the present con- 
text, such a construction will be achieved via the interac- 
tion of the system with specific classical radio-frequency 
fields, using the so-called Ramsey pulses Q. The den- 
sity matrix associated with |5'i2) is pi2 = |^i2)(5'i2| = 
|ciPpi + |c2pP2 + Api2. The crossed contribution: 
Api2 = ci C2|^'i)(^'2| + h.c. contains all the information 
needed to obtain the difference of the geometric phases 
acquired by the states and \^2) during an adiabatic 
quantum cycle. In some experimental schemes consid- 
ered in a forthcoming paper [H the geometric phase ac- 
quired by [^'2) will be a priori zero, and the measurement 
will yield directly the phase acquired by l^*!). 

A huge amount of work has been triggered by the pub- 
lication of M. Berry's paper [l[ on the phase acquired by a 
quantum state at the end of its adiabatic evolution along 
a closed cycle. Among the flurry of following papers, 
two selections have been presented in comprehensive re- 
views 6, 7] providing guides to the extensive literature on 
the subject. There also exist pedagogical presentations 
and a thorough review of manifestations giving many ex- 
aruples drawn from a great variety of physical systems 
[a-[lOl|. More recently. Berry's phase has become a topic 
of renewed interest, both theoretical and experimental, 
with regard to improved information processing and more 
specifically for its potential use in quantum computing 
[111 [l2| . At the same time, its unwanted manifestations 
in fundamental precision measurements, have to be care- 
fully investigated [l3l-[lH. Conversely, in other settings, 
Berry's phase might be just the right tool to detect still 
unobserved effects, such as parity violation in atomic hy- 
drogen p^.[l7j. 

This paper presents a detailed study of Berry's phases 
resulting from adiabatic quantum cycles performed by 
an arbitrary spin, interacting non-linearly with external 
electromagnetic fields. We assume that the Hamiltonian, 
governing the time evolution, involves both linear and 
quadratic couplings of the spin operator S. The addi- 
tion of a quadratic spin coupling generates new features 
which, as we shall see, can show up for S* > 1. The 
first purpose of this paper is to develop the formalism 
to calculate them, whatever the value of S, integer or 
half-integer. In addition, we shall give explicit numerical 
results for several spin values. Throughout this paper, we 
shall deal with the following set of quadratic spin Hamil- 
tonians: 



HiBit), E{t)) = 7s S • B(t) + 7Q (S • E(<))2 



(1) 



We have found in the literature only a few papers deal- 
ing with the Berrry's phase for a s pin submitted to a 
time- varying quadrupole interaction [isj , [loj . These deal 
with nuclear quadrupole resonance spectra in a magnetic 
resonance experiment involving sample rotation. How- 
ever, the absence of any magnetic interaction is assumed, 
which leads to level-degeneracy. These papers generalize 
the Berry's phase to the adiabatic transport of degener- 
ate states. In this case, the geometric phase is replaced 
by a unitary matrix given by the Wilson loop integral 
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of a SU{n) non-abelian gauge potential where n is the 
dimension of the eigenspace associated with a given de- 
generate quantum level |2C)I - [23 |. In contrast, in our work 
we are interested in situations where the energy levels of 
the instantaneous Hamiltonian stay widely separated. 

We impose on the model the extra constraint E-B = 0. 
As shown previously by one of us (C. B.) and G. Gib- 
bons 24, 25], the parameter space of i?(B,E) with 
E-B = is isomorphic to the space E(/9) relative to spin- 
one states, namely,the complex projective plane CP^. In 
[25} . the Berry's phase generated by the Hamiltonian ([1]) 
for S = 1, is found to be mathematically identical to the 
Aharonov-Anandan (A. A.) phase, if one uses an appro- 
priate parametrization of CP^. However, the physical 
contents are in general different, since, in contrast to the 
Berry's phase, the A. A. phase is not restricted to adia- 
batic quantum cycles. Hereafter, we derive an expression 
for Berry's phase relative to ([Ij for spins 5* > 1. In this 
case one loses the identity with the A. A. phase, which 
involves now closed circuits drawn upon the larger pro- 
jective complex spaces CP^^ . 

It is convenient to introduce the rotation R{t) which 
brings the x and z axes along the E and B fields 
and the associated unitary transformation U{R{t)). 
This allows us to rewrite the spin Hamiltonian ^ as 
i7(E,B) = 7sS U{R{t))H{\)UHR{t)), with the dimen- 
sionless Hamiltonian ■H(A) given by T-L{\) — -f 
\Sl?i~^. The parameter A = HjqE'^ /{'-/sB), combined 
with the Euler angles 6, ip, a relative to R{t) leads to a 
convenient set of coordinates for CP^. A quantum cy- 
cle within the time interval , T satisfies the boundary 
conditions e{Tc) = 6'(0), ^{Tc) = (p{0) + 2n<^7r, A(Tc) = 
A(0), a{Tc) — a{0) +naTr, where n^p and Ua are arbitrary 
integers. The eigenstates of H(A), '4>{\,m) are labeled 
with a magnetic number m by requiring that the analyt- 
ical continuation of V'CA, m) towards A = coincides with 
the angular momentum eigenstates \S, m). Thanks to the 
constraint E-B = 0, 'H(A) has several discrete symme- 
tries we shall use throughout this paper. The "m" par- 
ity (—1)'^""', associated with a 7r-rotation around i, is a 
good quantum number for 'H(A). Therefore 'H(A) can mix 
states ■i/'(A, mi) and ^{X, if and only if |toi — is 
an even integer. This selection rule implies that ip{X,m) 
can be written as a direct sum of even and odd blocks. 
This greatly simplifies the construction of the eigenstates 
of H(E, B), and the Berry-phase determination. 

After some manipulations involving Group Theory 
techniques. Berry's phase is written as a loop integral 
in the parameter space CP^: 

/3(m, X) — <j) p{m, A) (cos 9dip + da) — m{d(p + da) (2) 

where hp{m, X) is the average value of S along the B 
field direction. This quantity is obtained by taking the 
gradient of the eigenenergies of i?(E,B) with respect to 
jsB. We have performed an explicit calculation of the 
eigenvalues of 'H(A) for the spins S* = 2, 3 and 4 for val- 
ues of A running from to 2, 1.4 and 1.2 respectively 



(Fig.l). A superficial look at the above expression of 
/3(m, A) may give the impression that our result is, af- 
ter all, not so different from the case of a pure dipole 
coupling. However, there are specific features associated 
with the quadrupole coupling which appear more easily 
for special cycles where A and a are the only time varying 
parameters. The particular case 5* = 2, m = is espe- 
cially instructive to this respect. We recall that Berry's 
phase associated with a pure dipole coupling is vanish- 
ing for an arbitrary cycle with m = state. Performing 
the above simple cycle with the state S — 2, m — and 
assuming a constant quadrupole coupling, X(t) — 1, one 
obtains from a look at curves of Fig. 1 the quite remark- 
able result /3{m = 0, A) = /p(0, X)da — tt {mod 2tt). 

Our motivation for this paper is to make contact with 
experiment, having in mind the spectacular progress of 
atomic interferometry. One has then to face the prob- 
lem of the practical realization of a quadrupole coupling, 
having a magnitude comparable to the magnetic dipole 
one. For alkali atoms this is unrealistic if the Stark shift 
arises from a static E-field. The ac Stark shift [23| in- 
duced by a light beam is much more fiexible. However, 
in most cases, reaching values of A close to 1 requires 
the tuning of the beam frequency to be so close to an 
atomic line that the ac Stark effect induces an instability 
of the "dressed" atomic ground state. In typical cases 
this implies stringent constraints upon the duration of 
the Berry's cycle. As a result, the question of the valid- 
ity of the adiabatic approximation becomes crucial. This 
has led us to devote a full section of this paper to the 
precise evaluation of the non-adiabatic corrections. Our 
theoretical analysis is illustrated by numerical results for 
a few relevant cases. It provides guide lines for our forth- 
coming paper, devoted to experimental proposals. 

Our analysis of the non-adiabatic corrections proceeds 
in two steps. We first deal with the corrections associ- 
ated with the time derivatives of the Euler angles 9,ip,a. 
A convenient approach is the study of the Berry's cy- 
cle in the rotating frame attached to the B, and E 
fields, by performing upon the laboratory quantum state 
^{t) the unitary transformation = U~^{R{t))^{t). 
The corresponding Hamiltonian is obtained by adding to 
-/sBn{X) the extra term 758 - AB{t), where AB{t) is 
the magnetic field generated by the Coriolis effect. The 
longitudinal component AB// = — 7^^(cos0 (p + a)z, 
gives rise to a pure dynamical phase shift at the end 
of the Berry's cycle. As expected it incorporates /3(to) — 
m §{dip-\-da) at its lowest order contribution with respect 
to ?7 = —AB/z^-fsB). When a is the sole varying Euler 
parameter, the higher order terms in ry gives the com- 
plete set of non-adiabatic corrections. In addition, we 
show that the odd-order ones oc 7?^"+^ vanish exactly if 
one chooses for A the "magic" value A* (77) . On the other 
hand the transverse component ABj^, proportional to 
sm9 if, presents risks: it involves a linear combination of 
the spin operators and Sy and induces a mixing with 
opposite rn-parity states ( Am = ±1), possibly nearly de- 
generate unless stringent constraints are imposed upon A. 
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As regards the non-adiabatic corrections associated with 
A, we concentrate our attention to the ramping process of 
X{t) from A(0) = to X{T) ^ 1, the Euler angles keeping 
fixed values. Our method can be viewed as an extension 
of the rotating frame approach. We introduce the uni- 
tary transformation which makes H(A) diagonal within 
the 15*771) basis. The time evolution equation acquires 
an extra non-diagonal matrix A'H(t) oc A(t) which has 
the same effect as an rf pulse with sharp edges if X{t) 
increases linearly with t. This leads to rather large os- 
cillating non-adiabatic corrections exhibited in our work 
(see FigE]) . The standard procedure to tame them out is 
to give to A(t) a Blackman pulse shape [27 1. 

In the last section, we propose an "holonomic" proce- 
dure for the entanglement of TV non-correlated 1/2-spins 
(or N Qbits.) The basic tools are Berry's cycles generated 
by a Hamiltonian, formally identical to iJ(B(i), E(t)), 
except that, now, S is meant to be the total spin oper- 
ator of the N spins, X]t=i^j- "^^^ method is based on 
a known mathematical property: any non-correlated N 
1/2-spin state can be expanded into a sum of and 
Sz eigenstates. A given eigenvalue of will appear sev- 
eral times if iV > 2, but all the angular momentum states 
M have different symmetry properties under permuta- 
tions of the N spins, which leave H{'B{t), E(t)) invariant 
by construction. Thus the states ^5 , organized into an 
orthogonal basis, behave as if they were associated with 
isolated spins. An initial non-correlated state, written 
as the sum J2i s'^sm^sa/' transformed at the end 

of Berry's cycle into J2i s "-S m ^^P'''Pi^ ^)'^s m- With 
an appropriate choice of the cycle, we have been able to 
achieve maximum entanglement for iV = 4 , 5*2 = 1. 



II. THE INSTANTANEOUS EIGENFUNCTIONS 
OF /f(B,E) FOR A GIVEN ADIABATIC CYCLE 

In this section we construct the instantaneous eigen- 
functions of i/(B,E) for an arbitrary adiabatic cycle. 
The result is put under a form well adapted to the 
calculation of Berry's phase by group theoretical meth- 
ods. Our method applies to both integer and half-integer 
spins. 



in the frame attached to the fields B(t) and E(t), ignor- 
ing their time-dependence. For the explicit calculations 
to be performed in the cases S" > 1, it is convenient to 
introduce the dimensionless Hamiltonian 'H(A): 



H{B,E) = jsBhH{X); 

n{x) - Sz/h + x{s,/nf- 



(4) 



It is important to note that H(A) is invariant under 
three transformations. The first, 7i, corresponds to the 
reflexion with respect to the x, y plane, which changes 
Sx, Sy into ~Sx and —Sy, but has no effect on 5*2 

Ti Sx Ti ~ Sx , Ti Sy 71 = Sy , Ti Sz 7l = 5*2; . 

(5) 

The second transformation, 72, is the product of the re- 
flexion with respect to the ?/, z plane by the time reversal 
operation. The transformation of the spin operator un- 
der T2 obeys the relations: 



To ^ SxT2 



-Sx, T2 ^ SyT2 ~ Sy, T2 ^ SzT2 ~ Sz 



(6) 

The third transformation is a rotation of tt around the 
z axis, which changes Sx and Sy into —Sx and —Sy re- 
spectively, while leaving Sz unaltered. 

Tjj ^ SxTs = —Sx, TJ5 ^ SyTs = —Sy, 7^ ^SzTs — Sz- 

(7) 

(Note that 7i and % have the same effect on pseudo- 
vectors but opposite effects on vectors.) 

Let us now discuss some consequences of these invari- 
ance properties. To this end, let us introduce the eigen- 
vectors -0(771, A) of the Hamiltonian H{X), together with 
their eigenenergies £{m,X): 



HiX)iJj{m, A) = £{m, X)ilj{m, A). 



(8) 



All along this work the eigenenergies are supposed non- 
degenerate. In the limit A — )■ 0, the states ■il;{m, A) have 
to coincide with the angular momentum states \S,m): 
S'^\S,m) = h^S{S + l)\S,m) and Sz\S,m) = hm\S,m). 
The relations ([5]) imply that the quantum average of S = 
S/?i relative to \ip{rn, X ^ 0)), i.e. the polarization of the 
quantum state p(77i. A), lies along the z-axis: 



A. Instantaneous spin Hamitonian. Symmetry 
properties 

As a preliminary step, it is convenient to study the 
particular field configuration where E and B are along 
the X- and z-axes respectively, n = x, b = z, and write: 



H{B,E)^jsB Sz+7qE'S'x 



(3) 



where the term S^/3 = h'^S{S + l)/3, that plays no 
role in the calculation of Berry's phase has been omit- 
ted. In other words, H{B, E) is the spin Hamiltonian 



(0^(777, A)|S|V'(77l, A)) = {0,0, p(77l. A)}. 



(9) 



The invariance under 7i and T2 requires that the off- 
diagonal elements of the alignement tensor vanish: 

Aij = i(V3(777, A)|{i;,;,i;j}|Vi(777, A)) :^6ijAii . 

A third important consequence of the above invariance 
properties is that the eigenvectors ^(777, A) may be taken 
as real. This can be verified by noting that the matrix 
associated with 'H(A) in the angular momentum basis 
jS*, 777), using the standard phase convention, is real and 
symmetric, so that its eigenvectors are real. 
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Equation (O implies that the operator associated with 
the rotation of tt around z, T3 = exp(— itt S'z), com- 
mutes with the Hamihonian H{X). For convenience, 
let us introduce the "m-Parity " , V3 = exp(iS')73 ~ 
exp(i7r(S' — Sz))- The m-Parity of the angular momen- 
tum eigenstate state \S, m) is (— l)''^^™^. Within the an- 
gular momentum basis, the Hamiltonian 'H(A) can then 
be written as the direct sum of two matrices acting re- 
spectively on the states even and odd with respect to the 
operator 7^3: 

n{\)=Ueven{\)®Uodd{\)- (10) 

As a conclusion, we would like to stress that the field 
orthogonality condition plays an essential role in mak- 
ing the mathematical problem tractable for spins S > \. 
Otherwise, the problem would become rapidly compli- 
cated and any insight into Berry's phase physics gets 
blurred by the algebra. 

B. The instantaneous eigenfunctions of //(B(t), E(t)) 

Since we are going to use group theory arguments, it 
is appropriate to recall some basic facts about the rota- 
tion group in Quantum Mechanics. One introduces the 
unitary operator U (7?.(u, x)) associated with the rotation 
TZ{u, x) acting on the spin state vectors: 

U{n{u, x)) = exp , (11) 

where 72.(u,x) stands for the rotation R around the unit 
vector u by an angle x ■ This operator provides a unitary 
representation of the rotation group in the sense that it 
obeys the multiplication rule: U{Ri)U{R2) = [/(i?i-i?2), 
together with the unitarity relation U~^{R) = W(R). 
Applying the above rule to the case where i?2 is an in- 
finitesimal rotation, one derives the important relation; 

U-Hniu,x)) S C/(7^(^i,x)=7^(^i,x)•S, (12) 

which expresses the fact that the unitary transformation 
U{TZ{u,x)) rotates the spin observables. 

Let us associate with the orthogonal vectors E and B 
the trihedron (e, 6 A e, 6) which can be constructed by 
applying the rotation Re{0, 93, a) to the fixed coordinate 
trihedron {x^y,z)^ with 9,ip^a denoting the usual Eulcr 
angles: 

RE{e, V, a) = 7^(z, (^)7^(y, e)n{z, a). (13) 

To ensure the validity of the adiabatic approximation for 
the quantum cycles generated by iJ(B(t), E(t)), we shaU 
assume that the Euler angles are slowly varying func- 
tions of time. More precisely, we shall require that their 
time derivatives a^9,ip, - together with the time deriva- 
tive of the Stark- Zeeman coupling ratio A - are much 
smaller than the rate A Emin/fi, where AEmin stands 



for the minimum distance between the energy levels of 
H{B{t), E{t)). The adiabatic cycle within the time in- 
terval < t < T is specified by the boundary conditions 
involving the two finite integers and Ha '■ 

(y5(r) = (^(0) + 2n^7r; a(r) = a(0) + 7i„ tt. (14) 

Note that in contrast to the periodic variables ip{t) and 
a{t), X{t) and 6{t) recover their initial values at the end 
of the cycle. 

To proceed it is convenient to introduce the unitary 
operator U{R{t)) associated with the rotation R{t), given 
by: 

R{t) - 7^(z, ^{t))n{y, e{t))n{z, a[t)). (15) 

Since U{R{t)) belongs to a unitary representation of the 
rotation group, it can be written, using equation (1151) . as 
the following operator product: 

C/(i?(t)) = exp — • exp— • exp— i^a. (16) 
n h n 

Using the relation (|12l) with 7?,(u, x) = and the 

identity R{t)^^a • b = a • R{t)h, it is straightforward to 
derive the important relation: 

U{R{t))H{B{t), E{t))U^ {R) = HiBit), E(t)) (17) 

As a consequence, the wave functions ^(m, t) defined as 

*(to, t) = U{R{t))tij{m, A(t)), (18) 

are instantaneous eigenfunctions of i?(B(i), E(i)) with 
eigenvalues E(m, i3(i), i?(t)) : 

H(B{t),'E,{t))^{m,t) = U{R{t))H{B{t),E{t))'iJj{m,X{t)) 
= E{m,B{t),E{t))'^{m,t) (19) 
E{m,B{t),E{t)) = jsJiBit)£{m,X{t)). (20) 

We would like to stress that the instantaneous wave func- 
tions 5'(m,t) have, by construction, a well defined phase 
since, as shown previously, the state vectors i/'(m, A(i)) 
can all be taken as real. 



III. EXPLICIT CALCULATION OF THE 
EIGENFUNCTION PARAMETERS FOR 
BERRY'S PHASE EVALUATION 

We present for S — 2,3 and 4, explicit calculations 
of the eigenenergies £{m, A) and the polarization of the 
eigenstates, which, for symmetry reason, is along the B 
direction p(m. A) = {tlj{m, X)\Sz\^p{'rn, X)). (In this sec- 
tion, for the sake of simplicity we use a unit system 
where h ^ 1). We construct the matrix associated with 
H(S', A)within the angular momentum basis:|5m). It is 
convenient to rewrite the Hamiltonian in terms of the 
operators Sz, S'^ = Sx + iSy and — Sx — iSy 

1~L(^S, A) = Sz XSx , 

- Sz-^sj+j{{s+r + {s-r)+^s{s+i). 
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S=2 P(m,A) 




FIG. 1; Reduced energies £{m,X) and polarizations p(m, A) versus A for 5* = 2, 3 and 4. Since £{m,\ = 0) = p{m,X) — m, 
intersection of the curves with the vertical axis A = indicate their m labeling. The S = 2 curves exibit the 
symmetry relations £{m,\) — — f(— m, — A) and p(m, A) = — p(— m, — A), derived in the text. For m > there is a clear 
evidence for the degeneracies expected when A S> 1. 



Using textbook formulae, for the matrix elements of For S = I and 2 one finds easily: 
{S"^)^ in the \Sm) basis it is easy to write 'H(S', A) in 
matrix form. Using the invariance under the symmetry 
Ts introduced in subsection II. A, 7^(5', A) can be writ- 
ten as a direct sum of two matrices acting respectively 
upon the states \Sm) with even and odd values of S* — m, 
HiS, A) = 'HeveniS, A) ® HoddiS, A), One of order S and 
the other of order (5* + 1) depending on the parity of S. 



'HoddO',^) — 
^eDen(2,A) 



A/2 + 1 A/2 
A/2 A/2-1 

A + 2 ^/iJ2X 
^ij2\ 3A V372A 
VW2A A-2 
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rtod<il^,Aj - 1 5A/2-1 



(21) 



The Hamiltonians HeveniS, ^) and 'Hodd(3,A) as well as 
'Hei,en(4, A) and Hoddi'i, A) are given in Appendix B. 

The eigenvalues are obtained in a standard way by 
solving the two polynomial equations: 

det(a; l-He^en(2,A)) = - rjx^ X + AxX^ - Ax + 12X 
= 0, 

det(a; l-Hodd(2, A)) = x^ - 5xX + AX^ - 1 = 0. (22) 

The eigenvalue equations for S — 5 and 3 — 4 can be 
found in the appendix. Mathematica codes yield explicit 
expressions for the eigenenergies £{m,X). Although the 
formulae are rather complex, they lead to very accu- 
rate numerical results for all the values of A of inter- 
est. The accuracy is better than 10~^^. This result has 
been checked using numerical Mathematica codes, ob- 
tained directly from the matrix expression. To calculate 
the polarization p(m, A) we have used the fact that it is 
given by the derivative of the eigenenergy with respect to 
7s (Hellmann-Feynman theorem). Writing the eigen- 
value of H{B,E) as: E(m,B,i;) = jsB£{m,X), with 
A = ^^^^ , one gets immediately: 



p{m,X) 



dF.{m,B,E) 
h'ys d B 



£{171, A) — A 



d£{m, A) 
dX ■ 



(23) 



The results are given in FiglT] for both the energies and 
the polarizations. The eigenenergies £{m,X) can be la- 
belled unambiguously since for small values of A the 
eigenvalues are equal to m, up to corrections of the order 
of A^ . In FigH] there is evidence for near-degeneracies 
of opposite-parity level pairs: £{m, A) is approaching 
£{m — 1,A) when A ^ 1 for m = S — q, where q is 
an even integer such that < q < S . This was ex- 
pected since in this limit 'H(A) is dominated by the term 
oc S'^. The convergence is much slower for to < since 
the positive term A term has then to fight against a 
negative Zeeman effect. There are S degenerate doublets 
with £{m,X)/X ~ to^, as expected in the limit A ^ 1. 
In addition, by looking at the even-odd (or odd-even) 
pairs, one sees clearly that the pair £{S, X),£{S — 1,A) 
converges, without crossing, to the degenerate doublet 
£{S,X) ~ £{S - 1, A) ~ XS"^. The next lower pair ends 
as the degenerate doublet having the energy ~ A (5—1)^ 
and so on until one reaches the isolated level m — —S, 
which has no possibility other than converging towards 
the non-degenerate level with energy A x 0. While this 
behaviour is clearly exhibited in the simple case S = 2, 
only its two first steps are clearly apparent for S* = 3 and 
4 , but we have verified that the above picture is valid for 
all values of to by computing the ratios £(771, A)/ A when 
A> 1. 

A striking feature of the case S" = 2 in Fig[T] is the 
symmetry of the plotted curves under the transforma- 
tion TO — > —777. A —A, involving both the eigenen- 
ergies and the polarizations: £(m,X) — ^ —£{—m,—X), 



p{m. A) — > —p{—m, —A). To prove these symmetry prop- 
erties in the general case, it is convenient to perform 
upon the spin system a rotation of angle tt around the x 
axis, i?(a;,7r). By introducing the spin unitary operator 
t/(i?(i,7r)) — exp(— 7 5*2; tt), associated with the rotation 
R{x,tt), one can write the transformation law for the 
Hamiltonian H{X): 



U{R{x, Tr))n{X)U-\R{x, tt)) 



-U{~X). 



Applying U{R{x, tt)) to both sides of the eigenvalue equa- 
tion 'H{X)ip{m, A) ~ £{m, X)'ip{'m, A), one can write: 

U{R{x,T:))'H{X)i){m,X) = --H(-A)C/(i?(x, 7r))V3(TO, A) 
^ £{m,X)U{R{x,T:)^j{m,X). (24) 

The state U{R{x,TT)ilj{m, X) is an eigenstatc of H{—X) 
with the eigenenergy — ^(to. A). It coincides, up to a 
phase factor, with the eigenstatc 7/'(— to, — A), since the 
effect of the rotation R{x, tt) is to flip the spin component 
along the z axis. This completes the proof that: 

£{m, X) = —£{—m, —A). 

with the similar relation for p{m. A) as a consequence of 
Eq. 

In addition, the quantum-averaged polarizations 
p(to, a) satisfy the sum rule 



^p(to. A) = 0, 



(25) 



valid for any value of S and A. It is readily derived us- 
ing the fact that p{m. A) can be evaluated as the partial 
derivative of the energy with respect to 7s S (Eq. ([^ 1 
and noting the general structure of the polynomial equa- 
tion whose solutions provide the eigenenergies. 

To end this section we would like to say a few words 
concerning the state to = which is of particular inter- 
est for integer spins S >2. The fact that a second-order 
Stark effect can induce a polarization in an initially unpo- 
larized state was somewhat unexpected. To gain greater 
insight we have performed a perturbation computation 
in powers of A. Although one has to go to third order 
to find a non-zero effect the final result is given by the 
rather elegant formula: 

p(0. A) = ^X^S{S + 2) {S^ - 1) (1 + 0{X^)) . (26) 

It gives a rather accurate result for 5 = 2 if |A| < 0.4 but 
the domain of validity gets smaller for S=3: |A| < 0.12. 
Finally we would like to point out the smooth behaviour 
of p(0. A) for 5 = 2 in the vicinity of A ~ 1 and the 
remarkably simple exact values taken by the reduced en- 
ergy and the polarization at A = 1, namely f (0, 1) = 2 
and p(0, 1) — 1. At the time of writing, it is unclear 
whether this is a mere numerical accident or an indica- 
tion of something more profound. 
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IV. BERRY'S PHASE FOR ADIABATIC 
CYCLES GENERATED BY H{B{t),E{t)) ACTING 
ON ARBITRARY SPINS 

We start this section by a mini-review introducing 
the basic physical and mathematical features of Berry's 
phase concept. In particular we derive the general for- 
mula giving the Berry' phase in terms of the instan- 
taneous eigenfunctions of the time-dependent Hamilto- 
nian generating the adiabatic quantum cycles. Intro- 
ducing in this formula, the results of the previous sec- 
tion and relying upon Group Theory arguments, we per- 
form an explicit construction of Berry's phase relative to 
H(B{t),E{t)), valid for arbitrary spins. The final result 
is expressed as a loop integral in the CP^ space, using as 
coordinates the Euler angles and the parameter A. We 
show that for the case 5 = 2 a circular loop drawn upon 
a spherical subspace of CP^ lead to a loop integral very 
different from that of the case S = 1 which involves a 
magnetic monopolc Bohm-Aharonov phase. 



A. The Berry's phase as a physical observable and 
topological concept: a mini-review 

In this introductory subsection, we are going to follow, 
in several places, a presentation due to the late Leonard 
Schiff, in tribute of its memory. He gave in few pages of 
its venerable textbook [2^ a correct and precise treat- 
ment of the adiabatic approximation, involving a non- 
integrable phase. To study the adiabatic quantum cycles 
generated by the Schrodinger equation 



an{t): 



where H{t) ^ H{B{t), E{t)), 



it is convenient to expand the solution in terms of 
the eigenstates of H{t) : 

* (0 = W ^("' (27) 

n ^ ^ 

where x{n,t) = j{n,t) - [ dtiE{n, B{ti), E{ti)). (28) 

Jq 

The first term 'y{n,t) in (l28l) is a phase that vanishes 
for t = 0, but is otherwise arbitrary. Its value will be 
determined by the reasoning leading to the adiabatic ap- 
proximation. The second term, 

1 /■* 

<^D{n) = ^-J dhE{n,B{h),E{h)), (29) 

known as the "dynamical phase" , produces a contribu- 
tion which cancels H(t)^{t) in the wave equation. We 
shall take as the initial condition : $(0) = \E'(m, 0) or, 
in other words, a„(0) = Sn,m- In this case we can re- 
place the exact Schrodinger equation by the system of 
differential equations involving the expansion coefficients 



d 



ihdnit) = a„(t) {j{n,t) - ('^{n,t)\ih—'^{n,t))) 



Y,akit) exp «(x(fc, t) - x{n, i))(*(n, t)\i ?i^*(^, t)). (30) 



Since the adiabatic condition requires that a„(t) « 1 
whatever t, a necessary condition to ensure its validity 
is to cancel the coefficient of a„(t) in the r.h.s of equa- 
tion ((30)) . This is achieved if we make the following choice 
for the "gauge " j{n,t): 



* d 
dti(*(n,ti)|z?i— *(n,ti)). (31) 
ot 



Later on, we shall see that this non integrable gauge 
is a basic ingredient in the mathematical expression of 
Berry's phase in terms of the instantaneous wave func- 
tions \I'(m,t). The next step to validate the adiabatic 
approximation is to find appropriate conditions allow- 
ing the sum J2n^m Wn{t)\^ to remain below a predefined 
level for t > 0. This task will be performed in details 
in Sec. V but for the moment, let us assume that it is 
achieved. Within the adiabatic approximation, the solu- 
tion of the Schrodinger equation ih-^^(t) = H{t)^{t), 
with the initial condition $(0) = ^(m, 0), is then given 
by: 

^ADB{m,t) = 
t 



exp - \^j{m,t)~ J dtiE{m,B{ti),E{ti))j^{m,t).{32) 

We can now calculate the phase shift of the wave function 
^ADB{iTi,t) at the end of the adiabatic cycle: 



arg 



V «'az5s(to,0) 
1 



- 7(-,r) 



*(to,0) 



dtE{m,B{t),E{t)) , 



^ Jo 



dtE{m,B{t),E{t)) + I3{m), (33) 



where we have made Berry's phase /3{m) stand out on 
the r.h.s of the above equation. Using the equation pip , 
one gets immediately the basic formula giving /3(to) in 
terms of the instantaneous wave functions ^(m,t): 



I3{m) 



^ d 
dt{^{m,t)\i —'i!{m,t)) + 0(m), 



/o dt 
4>{m) = arg (^^(m, r)/^'(m, 0)) 



(34) 



It is crucial to note that the "dynamical phase " 0£)(m) 
and /3(m) obey different scaling laws under the trans- 
formation, B — ^i?, involving an arbitrary real pa- 
rameter ^, while keeping invariant the Euler angles and 
the dimensionless parameter A. If one remembers that 
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E(m, B{t), E{t)) = ^sB £{m, X{t)), one sees immediately 
that 4>£,{m) is multiplied by while I3(rn), being geo- 
metric, is unchanged. In principle, this scaling difference 
could be used to separate the two phases. However, a 
more practical way to isolate Berry's phase consists in 
measuring the phase for a given adiabatic cyclic evolution 
and that associated with the "image" circuit obtained by 
performing on the Hamiltonian parameters the transfor- 
mations: a{t) — — a(i) , Lp{t) — > —Lp{t), while keeping 
the other two unchanged. The two competing phases are 
transformed as 4)o(jn) — > (j3j){m) , (Urn) — >■ —/3{m), so 
that the dynamical phase can be eliminated by subtrac- 
tion. 

To end this mini-review, we would like to give, within 
the present context, a simplified description of the topo- 
logical interpretation of the Berry's Phase, due to Simon 
0. We have just shown that ^ADB{m,t) is a physical 
state obeying the Schrodinger equation within the adia- 
batic approximation. For our purpose, it is convenient 
to introduce the mathematical vector state <l>||(m,t) — 
exp (— i0£)(m, t))) $a_db(™, and to calculate the dif- 
ferential form ($11 ) = ($11 (m, i) | ^$|| {m, t)) dt taken 
along the adiabatic loop: 

d 

($|||rf$||) = {<^>ADB\-iw:<pD{m,t)^ADB)dt + 

" " at 

{^ADB\^H{t)^ADB)dt = 0. 

in 

The evolution of the state vector $|[ along the closed loop 
is then said to satisfy the "parallel transport" condition 
($l||d<i>l|) = 0. If the state is injected into the 

general formula for the Berry's phase, one immediately 
finds that 



/3(m) = arg( 



$ll(m,0)' 



(35) 



in the case of parallel transport. 

Let us now give a rather elementary introduction to 
the mathematical concepts behind the above notion of 
"parallel transport" applied to the evolution of quantum 
states. This arises rather naturally from a linear fiber 
bundle interpretation of Quantum Mechanics. The lin- 
ear fiber bundle associated with the quantum state space 
is constructed from the "base space" E{p) of the "pure " 
state density matrix p = \^){^\. Assuming for simplicity 
that the states 4* have unit norms, one finds that p sat- 
isfies the simple nonlinear relation = p. This implies 
clearly that E{p) has a non-trivial topology. The "fiber" 
is the one-dimensional space associated with a definite p. 
The vector states of the fiber are given by $ = exp(i0)^', 
where (p is an arbitrary phase and ^' a representative 
state of the fiber. The infinitesimal variation d$ during 
the quantum cycle is said to be "vertical" if it takes place 
along the fiber: = dr^Vi where dT is an infinites- 

imal c-number. Conversely, it is called "horizontal" or 
"parallel" if {<^\\\d<^\\) = 0. The fact that the Berry's 
phase can be viewed as a displacement respective to the 



fiber, associated with p(0) = p(T) - in our case a phase 
shift resulting from a parallel transport along a closed 
loop drawn upon the base space E{p) - emphasizes its 
topological character. 



B. The Berry's phase as a loop integral in the 
Hamiltonian H{'B{t),E{t)) parameter space 

Our starting point is the formula p4p giving Berry's 
phase for quantum adiabatic cycles associated with the 
instaneous wave function '^(m,t). 

The fundamental property of /3(m) is its invari- 
ance under the gauge transformation ^'(m, i) — > 
exp {j-f{t)) ^'(m, t), where f{t) is an arbitrary real func- 
tion. The density matrix of an isolated quantum system 
has clearly the same gauge invariance property. The adi- 
abatic approximation allows us to make a mapping of 
the Hamitonian parameter space onto the density ma- 
trix space. As a consequence, /3(m) could also be viewed 
as a line integral along a closed path drawn in the density 
matrix space. 

a. A group theoretical derivation To evaluate the 
expression p4p. it is convenient to write /3(to) in terms 
of ■(/'(m. A): 

/3(m) - 0(to) = ^ dt{^/j{m, X{t)\i ^VX™, ^(t))) + 

^ dt (^(m, \mU-\Rm ^ (C/(i?(i))) Him, A(t)).(36) 

Since [iplm, X{t)) is a real vector in the \S,m) basis, the 
first term of the integral can be written as: 

~{Kl-s [{s,pmm,\itm. 

If the cyclic condition A(T) — A(0) is satisfied, the line 
integral of this term is zero. If now one uses the explicit 
form of U{R{t)) in terms of the spin operator S (Eq. (|16l) . 
it is easily seen that the second term of Eq. p6p is a linear 
combination of time derivatives of the Euler angles: 

ihU-\R{t))-^^ {U{R{t)) + Dg + ipD^. (37) 

Let us consider Da- 

III 
Da = exp {-S^a{t)) ■ exp {-Sy9{t)) ■ exp {-S:,ip{t)) x 
n h h 

exp(-^5,^(t)) •exp(-j5y0(i)) -5^ •exp(-^5,a(i)). 

Since commutes with the exponential exp {—j-Sza(t)), 
Dr, reduces to: 



Da^U-\R{t))U{R{t))S,=S,. 



(38) 



For the two remaining terms such as the calculation 
is not so simple. Using the fact that 5*2 commutes with 
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exp {—j-Sz^it)), one can write 

= exp(-S'2a(t)) • exp(-S'y6'(i)) X 



•exp(--5j,0(i)) •exp(--5,a(i)). 



Using eauation (ll2l) . one can derive the commutation re- 
lation: 

S'^-exp {-^SyOit)) = exp {-^Sye{t))iSz cosO-Sa, sm9) . 

After repeating a similar operation to push 
exp{-j-Sza{t)) to the left of {S^cose - S^sin0), 
one arrives finally at the following expression: 

Dip — Sz cos 9 + sin 9{—Sx cos a + Sy sin a) . (39) 

In a similar way one obtains 

De — Sy cos a + Sx sin a, (40) 

but this term will not lead to any contribution to /3 be- 
cause of the boundary condition: 9{T) — 9{0). Since the 
quantum average of S relative to the state ip{m^ A) has a 
single component along z, we end up with the compact 
expression: 



/3(m) — 4>{m) 



p{m,X{t)){coa9 ip + a)dt . (41) 



We now have to calculate the phase shift (j){m) = 
arg {\I'(m, T)/^(m, 0)} appearing in Eq (p4)) . A pre- 
liminary step is to rewrite the unitary transformation 
U{R{t)) under a modified form: 



U{R{t)) = V{t) exp ( SzMt) + a{t)) 



V{t) = exp(-- Sz^it)] X 



exp\-^Sy9it)] exp(^Sz^{t)] .{42) 



The transformation law given by equation (1121) can be 
extended to any tensor operator SiSj 5*^ ... and conse- 
quently to any analytical function of Sy. Using the ex- 
tended law, one can write V{t) in the compact form: 



V(t) = exp \^--9{t){cosip{t) Sy + simpit) Sx) j . (43) 

Using the boundary conditions 9{T) = 9(0), X{T) = 
A(0),a(r) = a(0) + and ip{T) = Lp{0) + 2n^Tr, one 
obtains the relation: V{T) — V{0). Introducing for con- 
venience the notation U{z, u) = exp (— |- Sz u) it is then 
possible to rearrange 'i>{m,T) as follows: 

vI/(TO,r) = V{0)U{z,ip{T)+a{T))^im,m), 
= U(0) U{z, ip{0) + a{0)) U{z, Aip + Aa)?/'(TO, A(0)), 
with A(f + Aa^(p{T) + a{T)-ip{Q)-a{0), 

= (2n^ + na)n. (44) 



The next step is to prove that ^{m, A(0)) is an eigenstate 
of U{z., Aip + Aa) using its expansion in terms of the 
angular momentum state vectors \S,q), eigenstates of 5*2 
of eigenvalue hq: 

^(to,A(0))= CraAK0))\S,m-2n). (45) 

\m-2n\<S 



Using equations (j44|) . one sees immediately that, as a 
consequence of the quantum cycle boundary conditions, 
each state \S,m — 2n) appearing in the above sum is 
an eigenstate of U{z, Aip + Aa) with the eigenvalue 
exp {—imTT(2nip + Ua)). This leads to the basic relation 
connecting 'I'(r7i,r) and 'l'(m,0): 

5'(m,r) = exp (— i rn 7r(2 n,^ -|- tIq)) X 

U(0)C/(z,^(0)-H«(0))7/;(m,A(0)) 
exp (-im (A(/7 + Aa)) ^'(m, 0). (46) 

One immediately obtains the phase shift 0(m) appearing 
on the r.h.s of equation 



(47) 



By combining the equations (|4ip and (|47p. we arrive 
at the expression of Berry's phase for an arbitrary 
quantum adiabatic cycle generated by the Hamiltonian 
H(B{t),'E{t)) of equation (1) with E{t) ■ B{t) = 



0(m) = arg{*(m,T)/\E'(TO,0)} 

= -m{^{T) ~ ^{0) + a{T) - a{Q)) 

^ -m (if{t) + a{t)) dt . 



f3{m) 



[m — p(m, \{t)) cos 9{t)] (p{t) dt 



[m - p{m, X{t))] a{t) dt . (48) 



b. The Berry's phase as a loop integral of an Abelian 
Gauge Field The Berry's phase /3{m) given by equa- 
tion (j48p can be written as a loop integral around a 
closed curve C drawn in the parameter space CP^, using 
A, 9, ip, a as coordinates: 



j3{m) — ~ (p ((m — p(m. A) cos 9)dip + (m — p{m^ X))da) . 

(49) 

The above formula suggests that /3(to) could be written 
as a loop line integral of an Abelian gauge field Ai{xj), 
defined in the parameter space. The non-vanishing com- 
ponents of the gauge field candidate, A^p and Aa, are 
then given by: 



A,, 



-m + p{m,X) cos 9, Aa — ~m + p{m,\). (50) 



Thus the Berry's phase looks like the Bohm-Aharonov 
phase |29j, for the Abelian gauge field {A^p.Aa)). 



P{m) 



<j> {A^ dip + Aada). 



(51) 
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FIG. 2: Representation of the gauge field component Aci{9) 
versus the associated 9 parameter of the non-hnear Hamilto- 
nian in the two cases 5=1 (red curve) and S — 2 (black 
curve). The difference between the two graphs illustrates the 
fact that for S > 1 there is no longer a one to one corre- 
spondence between the Hamiltonian parameter space and the 
density matrix space upon which is drawn the quantum cycle. 
The non-linear spin Hamiltonian generates a new geometry 
revealed only by spins larger than one. 



ponent of the gauge field Aa for two different cases. For 
S' = 1 , m = 1 , the gauge field is given by the same expres- 
sion as tfie gauge field associated with a linear Hamilto- 
nian Aa{9) = COS 9 — 1. However, the situation is totally 
different for the case 5 = 2, m = 0, where Aa{9) exhibits 
the peculiar shape shown in Fig. 2. This implies that the 
geometry involved in Berry's phase for the sub-manifold 
differs from its usual interpretation in terms of the solid 
angle defined by a closed loop drawn upon a 2D sphere. 
The origin of this phenomenon lies in the fact that the 
density matrix space for S* = 2 is not CP^ but rather 
CP"*. Unlike the case 5 = 1, there is no longer a one 
to one correspondence between the parameter space and 
the density matrix space. It is thus not surprising that 
the geometry involved in Berry's phase for a non-linear 
Hamiltonian should differ from the linear case, even for 
the simple cycles involving only the variations of the two 
parameters A and a. 

V. NON-ADIABATIC CORRECTIONS IN 
BERRY'S CYCLES USING THE ROTATING 
FRAME HAMILTONIAN 



What remains to be proved is that the above expression 
is indeed invariant under the gauge transformation of the 
Abelian field: 

d d 

A^ A^ + —g{Lp,e,a,\), A^ Aa + -^g{'f,9,a,X). 

The crucial point is that g(ip,9,a, X) must be a single- 
valued function of B and E E, like the gauge field 
itself. This implies that g{(p, 9, a. A) must be a peri- 
odic function of ip and a, with periods respectively of 
27r and tt. Remembering the boundary conditions (jl4p 
for the quantum cycle C , one finds immediately that the 
gauge contributions to (3{m), g{(p{0) + 2n^Tr,9{0),a{0) + 
n„7r,A(0)) - g(v3(0),6'(0),a(0), A(0)), do indeed vanish. 
This completes the identification of /?(m) as a Bohm- 
Aharonov phase. 

c. Berry's phase geometry generated by non-linear 
spin Hamiltonians for S > 1 For the sake of simplicity, 
in the next sections we shall concentrate on Berry's cycles 
where A and a are the only time- varying parameters. The 
question then arises as to whether we shall not lose in this 
way most of the new features introduced with the non- 
linearity of H{t). In order to make clear the cause for our 
concern, let us perform the change of variables upon the 
two left-over parameters which map the associated 2D 
manifold onto a 2D sphere: A = —2 cot{9) with < 9 < tt 
and tp ~ 2 a. In ref. (23 | we have given an explicit form 
of the CP^ metric obtained with this kind of coordinates. 
For the sub-manifold associated with (p = 9 = it turns 
out that this metric, ds^ — j ^dt^^sin^ 9 -\- d9^^ is iden- 
tical to the one associated with a 2D sphere having a 
radius 1/2 (this factor results from the choice: if = 2 a.) 

To answer the question, we can now consider the com- 



For quantum cycles of finite duration T, the time 
derivatives of the Hamiltonian parameters cannot take 
arbitrary small values. For instance, in the experiments 
discussed in a separate paper [SJ , the quadratic spin cou- 
pling generated by an ac Stark effect induces an atomic 
instability when Awl. As a result, the measurement 
of Berry's phases studied by interferometry experiments 
on alkali atoms will require a good control upon the non- 
adiabatic corrections arising from the time dependent ex- 
ternal B, E fields. 

We shall consider separately the effect coming from 
the periodic parameters a and ip in subsections V.A and 
V.B, and the non-periodic one A in subsection V.C. Our 
method involves the time-dependent Schrodinger equa- 
tion in the rotating frame. We shall use an adiabatic 
approximation within this frame, the slowly varying pa- 
rameters being, then, the time derivatives of the periodic 
parameters a and ip. 

We set = U{R{t))^{t), and we write the wave 

equation relative to <f>(t) in the rotating frame: 

(^,nj^_H(^B{t),E{t))^ 

-U-\R^t))^hJ^U{R{t))m = 

-{D^ a + D^ip + De9) $(t) = 7s S • AB{t) $(t). (53) 

The rotating frame Hamiltonian can be written as 

H{t) = H{B{t),E{t)) + 7sS • AB{t), (54) 

where AB{t) is an additional effective magnetic field gen- 
erated by the Coriolis effect. It can be decomposed into 
a longitudinal component AB / / (t) and a transverse one 
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A'B±(t). Using the explicit expressions of Da and D^p 
in subsection III.B and assuming here -for the sake of 
simpHcity- that 9 — 0, one arrives at the following ex- 
pressions for of the effective field: 

AB// = — J {cos 9 ip + a) z 

AB_L = — J sin 9 ip {— cos a X + sin a ij) 

= sin9 ipTl{z,—a) X. (55) 

The role of those two components are going to be exam- 
ined separately. 



A. The Berry's phase and its non-adiabatic 
corrections involving the longitudinal effective 
magnetic field {B{t) + AB/ /{t)) z 

We consider, first, the part of the Hamiltonian govern- 
ing the evolution of $(/:) associated with AB//, which is 
the sole present if = 0. This part can be written very 
simply in terms of H{B,E) 



H//{t)=H{B{t)+AB,/{t),E{t)). 



(56) 



To proceed it is convenient to introduce the dimensionless 
parameter ry : 



AB 



' 1 1 cos 9 ip + a 



(57) 



B jsB 
Using B = {l — ri)B, the eigenvalues of H//{t) are simply 
E{m,B + AB//,E) ^ -isBh£{m,\) 



= Is Bh{l - r])£{m, 



X 



1-77 



1(58) 



Since the quantum cycles generated by H/ /{t) are devoid 
of any topology, we can expect Berry's phase relative to 
the laboratory frame to be buried inside the dynamical 
phase of H{t): 

4>D//{X,v,m) = - [ dt-isB{l-ri)£{m,-^). (59) 

Indeed, it is easily seen that Berry's phase is given by 
the first-order contribution in the 77-series expansion of 
the above integral : 



dt p{m, X) (cos 9(p + a), 



/3(rn) — (pim), 



(60) 



where we have used equation (l23l) of sec. III. The phase 
(j){m) is recovered by returning to the laboratory frame. 
The even-order contributions cx rj^^ are eliminated by 



subtracting term by term those coming from image cir- 
cuits associated to opposite signs of the 77-parameter (i. 
e. — >■ —if and a —a). 

We now show that it is possible to find a magic value 
of A = A* (77) that leads to a cancellation of all the contri- 
butions of the dynamical phase cx 77^^+^. More precisely 
let us define 

A/3// = i (4>D//{X,r],m) - 4>d //{X, -ri,m)^-l3{m)+(f){m) 



dt{cos 9 ip + a) Ap(rn, A, 77), 



(61) 



where Ap(777,, A, 77) ~ 

i-(^77+l)f(7n,- 



V 



) - (77 ^- -7/) -p{m,X). (62) 



The above formula seems to suggest that the inte- 
grand in the expression giving A/3//{m) is singular when 
7] — ±1. In fact it is easily seen that this is not the 
case. Indeed, the ratios £(m,X)/X converge towards the 
eigenvalues of (Sx/h)"^ : 0,1, 4. ..5*^ , as A — )■ 00. As a 
preliminary step, let us consider the lowest-order expan- 
sion of the r.h.s of Ap(m, A, 77) with respect to 77 



Apf/{m) = 



dt {cos 9 ip+ a) q{m,X)rf' , 



q{m,X) = —i^X^£"'{m,X) + 3X^£"{m,X)j{63) 

Fig|3]shows the variations of q{m, A) versus A in the cases 
5 = 2 and 4, for tti = 0. A quite remarkable effect, clearly 
visible on Figl^l is the vanishing of q{0, A) for a particular 
value of A , in the vicinity of the maximun of -^p{0, A). 
This should allow for rather accurate determinations of 
/3(0), as explained in the caption of FiglH 

In the case where a is the sole time-dependent parame- 
ter we have found, that a similar property holds in fact to 
all orders in rj. More precisely, we have shown that there 
exists for a given value of 77 a "magic" value X*{S, rj) such 
that the non-adiabatic correction to p{m, A) cancels ex- 
actly. 

Ap(0,A'^(5,77),77)= 0. (64) 

The magic values can be accurately represented by the 
polynomial fits given below: 

A* (2, 77) = 0.838213- 0.083782377^- 0.0431478774- 

0.02318877/6 - 0.020798677^ 
A* (4, 77) = 0.509982- 0.0900927772- 0.034998577'*- 

0.04364957/6 + 0.0373634?/^ (65) 

As it is apparent upon the above formulas, the "magic" 
values X*{S,ri) are slowly varying functions of 77 within 
the interval < 77 < 0.5. We have verified that they 
both satisfy the equation with a precision better 
than 3 x 10"^. In conlusion, the above results open the 
road to measurements of the Berry's Phase (3{0) for S = 2 
and 5 = 4, for quantum cycles where a is the sole varying 
Euler angle, under conditions where the non-adiabatic 
corrections can be kept below the one ppm level. 
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FIG. 3: Comparison of Berry's phase to its lowest-order non-adiabatic corrections, when a is the only time-dependent Euler 
angle, for a quantum cyle involving the spin configurations 5 = 2, or 4 and m — 0. The Berry's phase is given by /3(0) = 
dtp{0, X{t))a{t). Using the results of subsection V.A, one finds that the non-adiabatic correction is given by : A^'^'/3(0) — 

dt g(0, \{t))a{t) ( -y^B(t) T^^s explicit expression of g(m. A) in terms of the second- and third-order derivatives of £{m, A) 
with respect to A can be found in Eq. (|63p . One sees clearly that, for 5* = 2, if the parameter A is required to stay within the 
small interval 0.8 ± 0.2, then A^^'/3(0) //3(0) will be < 10~^ provided a{t) is kept smaller than one tenth of the Larmor angular 

frequency 'ys B{t). Similarly, for 5 = 4, there is a value of A, for which q{0, X) cancels. The partial derivative ^p(0, A) is 
governing the response of /3(0, A), say, to a small variation of the B field, are nearly at their maximum when q{0, A) is close to 
0. Note also the fourfold increase of the sensitivity for 5 = 4 compared to S = 2. 



B. Non-adiabatic corrections induced by the 
effective transverse magnetic field ABx(t) 



This subsection is divided into two parts. In the first 
we calculate the corrections to the Berry's phase of the 
laboratory frame that are odd under ?7-reversal. These 
arise from the modification of the dynamical phase in- 
duced by ABj^(t). Then we turn to the small Berry's 
phase in the rotating frame associated with the non- 
trivial topology of the quantum cycles. It arises from the 
time variation of the effective field B{t) z + AB{t) acting 
in the rotating frame. Let us stress again an important 
feature of the transverse contribution, 75'S • ABj^(i) to 
H{t). Being proportional to cosaSx — sinaS'j,, it mixes 
states of opposite m-parity which can be nearly degener- 
ate, unless severe restrictions are imposed upon the do- 
main of variation of A. These constaints are directly read 
off from Fig.l. For S = A,m = 1,2 and 3, only negative 
values of A are allowed, —2 < A < 0. For S* = 2, to = a 
larger interval can be used: |A| < 1. 



1. Correction to the dynamical phase 

It is convenient to introduce the dimensionless param- 
eters /X and Jl and the dimensionless operator S = ?i~^S: 



AB± sin dip _ /i 
= 7^ = TT' Ai = 



(66) 



B jsB 

We rewrite H{t) as follows: 

H{t) = -fsBh{n{X) - cosaJ:^ + smaJ:y)). (67) 

To proceed, it is convenient to introduce the perturbation 
expansion of the eigenenergies of H{t): 

ECX,tl)^lsBhi£{m,~X) + il^£^f\m,\)+Oifl^)). (68) 

The second order energy shift caused by the Jl- 
contribution to H{t) is given, to all orders in rj, by: 



(to,A)=^ 



|(V'(n, A) I cosa Sa; — sma'Ey\^p{m, A))p 



\n—m\odd 



£{m, A) — £{n, A) 



(69) 

The symmetry properties of H{X) (Sec. II) have two con- 
sequences i) only \n — to | -odd terms contribute to the 
sum and ii) the cross terms involving E^jEj^ cancel out. 
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This means that there is no contribution proportional to 
sin a cos a. We are left with the terms prop to cos and 
sina^. We shall assume that the velocity d is a slowly 
varying function of t during the cycle, so that cos and 
sinck^ can be replaced by their average values of 1/2. 

At this point we have found convenient to solve nu- 
merically two auxiliary problems, namely the search of 
the eigenstates of the two Hamiltonians: 

Hi (A, /i) = 'H(A) — jiYji for i — x and i — y. 

Let £x,y{X,fJ') be the individual eigenenergies of 'Hx,y 
which are even functions of /i. To lowest-order in /i, they 
take the form: 



rotates about the z axis with the angular velocity a. As 
above we shall limit ourselves to the second order con- 
tribution with respect to /i. We proceed by writing the 
eigenfunctions of 'H(i) in first order expansion with re- 
spect to the parameter pt: 



i^im, A) - /i?/>(i'(m, A) -I- 0{fl^). 



(74) 



In a way similar to what we did before, it is convenient to 
introduce the first order expansion of the eigenfunctions 
Tl'x,yim,X) of 'H^,j;(A, /x). 



£x,y{X,^l) = £{X) + E^^hX) + 0(/i^), 



(70) 



^x,y{m,A]~ 2^ £(m,A)- 



\n—7n\odd 



{'ip{n,X)\Y.x, y\'ip{m,X)) 
£{n, A) 



7/>(n, A). 



where £x,y{X) are obtained by an interpolation towards 
^ = 0. The overall second order energy shift of Ea. ([M)) 
can be expressed as: 



By using Eq. (j67|) . one finds that the first-order contri- 
bution to the H{t) eigenfunction (divided by jl) can be 
written as: 



:f (™,A) = i(42)(A)+42)(^))^ (7^) ^W(m,A) = -cosa^W(m,A) + sina^W(m,A). (75) 



We have now all we need to write the contribution to the 
dynamical phase associated with H{t) to 2nd-order in /i 
and all orders in rj: 

dt-isBjj?£f\m,X), 



/,(2) _ 

■'d± — 



dt-fsB 



to 1-r] - ■ L-ri' 

In practice, we shall concentrate on the contribution to 
the dynamical phase relative to H{t) of the order of fi^ rj, 
which is easily derived from the above equation. In order 
to highlight the close connection with Berry's phase we 
replace r] in the integrand by its explicit expression of 
Eq.dSll): 



f\m,^). (72) 



A 



'(/3-|-d)p(^)(m,A), 



with p^^^m,X) = {l + X^)£^j'\m,X)). 

OA 



(73) 



Despite the sign difference in the A derivative, our nota- 
tion underlines the analogy of the above correction with 
the Berry's phase, but the presence of the factor jj? oc 
spoils its geometric character. The quantity ^'•^■'(to. A) 
has been plotted in Figure |4] for the cases 5 = 2, m = 
and m — —1 over the interval 0.7 < A < 1.2. As 



expected, the non-adiabatic corrections are 



< 



for 



m — 0,-1. We have performed a similar computation 
for m — 1 and found, in constract, that p^^^m, = 1,A) 
blows up to values ^ 100. 



2. The Berry's phase associated with AB± 

The full Hamiltonian H{t) /f//(i) +758 • AB^, gen- 
erates quantum cycles endowed with a non-trivial geom- 
etry since the effective transverse magnetic field ABj^ 



It is then easily seen that Berry's phase appears only to 
second order in and involves the familiar time derivative 
product: dt{^p'^^^m,X)\ih§-^'tP<^^\m,X)), which IS pro- 
portional to d. If the other parameters in the integrand 
vary slowly during the closed cycle, the only terms sur- 
viving are those proportional to cos^ a and sin^ a which 
can be replaced by their average values of 1 /2. This leads 
to the following correction to the Berry's phase 

A(2)/3j.(A,/i) - r dtp'aQ{{i:^y'\m,XM'\m,X))}, 
Jo 

(76) 

a result valid to all orders in r/. It is important to note 
that, unlike the case of the dynamical phase (Eq. 72), 
this correction invoves the product of fi^ by d. There- 
fore, it is the contribution odd in 77 that is eliminated by 
the parameter reversal, so that the dominant contribu- 
tion corresponds to the limit ?7-^>0 i.e A— >A,/i— >/i. 
The quantity Cx,y{m, X) — 3{(Vi^"^''(m, X)\ij^x\m, A))} is 
plotted together with p^'^\m,X) in Fig. 5, leading to the 
same conclusions as those reached at the end of the pre- 
vious paragraph. 



C. Non-adiabatic corrections associated with 
ramping-up the quadratic interaction 

In the previous subsections Berry's phases were in- 
duced by cyclic variations of the periodic Euler angles 
performed in an initial state V'(*S'rn; Aq). In fact to pre- 
pare this eigenstate one has to start from an eigenstate 
of S and Sz at a certain time t=0. During the time 
interval < t < T, the evolution of the spin system 
is governed by 'H(A(t)) = S^/h + iX{t){Sx/hf. Writing 
X{t) — Ao g{t/T), we assume that g{s) rises "slowly" from 
to 1. We now ask asks whether, for a given values of 
the rising time T, there might be a better choice for g{s) 
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S= 2 m=0 




S=2 m=-l 



-0.30 



-0.35 



1.0 1.1 1.2 




FIG. 4: Two non-adiabatic corrections to Berry's phase associated with the effective transverse magnetic field ABj_(t) acting 
in the rotating frame. The first, <l>'^^ = dt {cosOip + a)p^^\ni, X), is the modification of the dynamical phase of H{t) 
associated with AB_L(t). The second, A*^'/3x (A, /i) = dt il^aCx,y{m, X), is Berry's phase in the rotating frame induced by 
the rotation of the transverse field, = ^b* ' ~ T^)- Note that both corrections are not geometric phases because of the 
presence of the factor fi^. The results for the cases 5 = 2, m = and -1 are of the same order of magnitude. They can be used 
to make Berry's phase determination precise at the order of fi*. 




FIG. 5: Effect of taming the non-adiabatic oscillations generated by a too abrupt application of the quadratic spin interaction, 
thanks to the use of a Blackman pulse-shaped time-derivative, oscillation taming (OT) procedure. The diagram shows non- 
adiabatic corrections to (Sz) at the end of the transition from A(0) =0 to A(T) = Ao, starting from an initial state which is an 
angular momentum eigenstate 15* = 2,m) (left, Ao = 1, m = —1; right, Ao — A*(0) = 0.838, m = 0). The big black dot stands 
for the adiabatic prediction for (Sz)- The red curves are for a linear increase of A(t), while the black ones correspond to the 
OT choice, i.e. X{t) = /(f )//(0), where /(s) is a Blackman pulse. The rather strong oscillating non-adiabatic corrections to 
{Sz} are clearly swept away by this latter choice, provided BT > 25. 



16 



than just g{s) = s when Aq ~ 1, if one wishes to keep the 
non-adiabatic corrections below a predetermined level . 
In the present subsection, we have chosen 1/(75 i?) as 
the time unit: all the time dependent physical quantities 
and their time derivatives are functions of t — t^s B. 

1. The two-level spin system S — 2,m — ±1 

To get some insight into this problem, let us first 
consider the simple case of the adiabatic evolution of 
the state S = 2, to = ±1 governed by the Hamilto- 
nian H2(i) = 'Hodd{'2,X{t)) defined by the 2 x 2 ma- 
trix of Eq. (l2T]) . By performing the change of variables 
C = arctan(|A) with — 7r/2 < C < 7r/2, one can write 
'H2{t) under the convenient form: 

5 

■H2 (t) = sec ({ax sin C + cos C) + ^ tan(C)cro , (77) 

where ctq is the 2x2 unit matrix and cr^^y.z are the 
standard Pauli matrices. We define $(2to; t) as the 
state vector obeying the Schrodinger equation relative 
to 'H^it) and satisfying the initial condition $(2to;0) = 
|2, to). We introduce the unitary rotation matrix V2(C) — 
exp — i^(Tj,/2. This is the rotating frame Hamilto- 
nian 7i2{t) which governs the evolution of $(2m;i) = 
V2~^(C) <&(2TO;t). After some simple algebra we arrive at 
the familiar expressions: 

5 1 ■ 

= 3 tan(C)(To + sec{C)a, - -(oy. (78) 

If C is assumed to grow linearly during the time interval 
< r < ^sBT, 'H2{t) can be identified with the Hamil- 
tonian describing a square Ramsey pulse. The same anal- 
ysis can be repeated for the 5* = 1,to = ±1. The result 
looks very similar : the spin dependent part of the ro- 
tating frame Hamiltonian 'Hi{t) is formally identical, but 
with one significant difference, namely tan^ = /2. 

In section VI, we shall study the time evolution of a 
four-spin system governed by the direct sum of the two 
Hamiltonians: 'Hi(t) © 'H2{t) relative to the same value 
of A. So we will have to use the two different expressions 
of C, for S=2, C2 ^ 6A/(9A2 + 4), and for S=l, C'l -> 
2A/(A^ +4). If \[t] is rising linearly, the C,\t] pulses are 
replaced by trapezoidal Ramsey pulses with sharp edges 
leading also to non-adiabatic oscillating corrections. 

There is a standard way to wash them out which is 
used in nuclear magnetic resonance and atomic interfer- 
ometry experiments. This consists in using the so-called 
Blackman pulse shape [l^]. Adapted to the present con- 
text, it leads the following choice of \{t): 

\BM{t) = Xof{l)/f{0), 

where f{s) = 0.42 - 0.5 cos (27rs) -I- 0.08 cos (47rs).(79) 



We have solved numerically the Shrodinger equation for 
Ao = 1 for the two different choices, Blackman and trape- 
zoidal pulses for X{t). The important quantity for our 
purpose is (5*2 (T)), since it governs the magnitude of 
Berry's phase induced by the cyclic variation of the Euler- 
angles. It is important to note that (^2) is affected to 
first-order in X(t) in contrast with the dynamical phase 
which is modified only to second order. 

We have plotted {Sz{T)) in Fig. [5] (left-hand curves), 
as function of ^sBT- The big (black) dot on the or- 
dinate axis shows the adiabatic prediction. The use of 
the Blackman pulse for X{t) leads to spectacular conver- 
gence towards the adiabatic limit: for jsBT > 25 the 
non-adiabatic correction to (5*2 (T)) plunges down below 
the 0.02 % level. We have also computed the phase shift, 
which is purely dynamical in the present case. As ex- 
pected, the non-adiabatic corrections are much less sen- 
sitive to the time dependence of X{t) and we find that 
in both cases \(j)oiexact) — (j)oiadiab)\ < 0.008 when 
'YsBT > 25. Similar calculations for S = \m\ = 1 and the 
initial value to = —1, for Aq = 1 indicate that corrections 
(not shown) are an order of magnitude smaller. 

2. More than two levels 

We now extend the foregoing "rotating" frame method 
to non-adiabatic processes involving the mixing of more 
than two levels. As noted in section II. A, the Hamilto- 
nian niX(t)) = S,Jh + X{t){Sx/nf is described in the 
angular momentum basis |5to) by a real symmetric ma- 
trix. This implies that its eigenstates ip{Sm] A) can be 
written as real vectors. An evident consequence is the 
following identity, obtained by taking the derivative of 
the normalization condition of the real eigenvectors: 

(^(5to;A)|^^(5to;A)) =0. (80) 

It implies that the adiabatic phase is purely dynamical. 
Let us introduce now the transformation matrix: 

T/5(A)=^|^(5to;A))(5to|. (81) 

One can directly verify that Vs{X) transforms H{Sm; A) 
into a diagonal matrix within the basis |S'to), by writing 
Vs{X)\Smi) = J2m'4'(Sm;X){Sm\Smi) = '0(S'toi;A). 
The time evolution in the "rotating" frame is then gov- 
erned by the following Hamiltonian: 

iisit) = V^HiSm;X)Vs-iXV^^Vs, 

= ns{x) + Ans{t). (82) 

The expression HsiX) = J2m^i^''''^T ■^)\^^) {Sm\ fol- 
lows from eq. ([5T|) . while the identity implies that 
AHsit) is non-diagonal: 

Ansit) = J2 l^"^2) (5toi| X 
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By using an identity obtained from the time derivative 
of the eigenvalue equation ([8|), we can rewrite the r.h.s of 
the above equation as: 



AHsit) - -iX{t) X J2 



\Sm2) {Smi 



£{Smi; A) — £{Sm2; A) 



mi 

(^(5mi; X)\{S^/hf\i^{Sm2; A)) . (83) 

The above formula can be considered as an adaptation 
of a more general one to be found in |28ij and refer- 
ences therein. One must stress that energy denominators 
£{Smi] A) — £{Sm2] A) involve levels for which toi— 7712 is 
even integer. For instance, in the cases 5 = 2, m = 0, ±2 
the two energy differences involved are varying slowly 
within the interval < A^ Ao = 0.838. It follows that 
the time dependence of A7is{t) is dominated by —i X{t) . 
Let us assume that X{t) is increasing linearly within the 
time interval < t < T : X{t) = Xot/T. The effect of 
AUsit): as in the previous cases, is equivalent to a r.f. 
trapezoidal pulse with sharp edges. In the right-hand 
curve of Fig. 5, we have plotted the average polarization 
as a function of the rising time T for X(t), towards the 
value Ao = A*(0) = 0.838. There is a a very rapid damp- 
ing of the oscillating non-adabatic corrections, passing 
below the level of 0.2% for 7s T > 25 when X{t) has a 
Blackman pulse shape. This is in strong contrast with 
the case of a linearly rising X{t) where, for 75 BT > 25, 
the non adiabatic oscillating corrections have an ampli- 
tude of about 10%, decreasing only very slowly for longer 
rising times. The conspicuous chaotic behaviour in the 
3-level case reflects the fact that the two frequencies in- 
volved have in general an irrational ratio. 



D. Concluding remarks 

Let us sum up the results of this last section. We have 
shown that if X(t) is a linear function within the finite 
time interval < t < T, vanishing elsewhere, oscillat- 
ing non-adiabatic corrections are generated by the sharp 
jumps of X{t) at t = and t — T. We have found that 
the remedy is to choose for X{t) a Blackman pulse shape. 
It is clear that one should use the same prescription for 
a{t),ip{t)^ for taming analogous unwanted oscillations. 
By solving exactly the Schrodinger equation governed by 
the Hamiltonian H//{t) of Eq. ([55]) . as illustrated in Sec. 
VI and in a separate paper [5] on two definite examples, 
we have indeed found that this remedy works beautifully. 

It should be stressed that the near-degeneracies of 
opposite-parity level pairs {£{m, X),£{m — 1, A)), which 
occurs when m > and A > 0, do not affect the validity 
of the adiabatic approximation for Berry's cycles if a(t) 
and X(t) are the only time-dependent parameters. This 
simplification is a consequence of the symmetry proper- 
ties of the Hamiltonian H{'B{t), E(i)) resulting from the 
E and B orthogonality. 



Furthermore, choosing for X{t) the "magic" expres- 
sion X*{-^^) given by Ea. (|55|) makes the difference of 
the Berry phases relative to the mirror cycles (a(t), X{t)) 
and (— a(t),A(i)) free of non-adiabatic corrections: all 
the correcting terms oc {a{t)/^s B)'^'^^^ vanish whatever 
g > 0. 



VI. HOLONOMIC ENTANGLEMENT OF N 
NON-CORRELATED 1/2 SPINS 

A. Introduction 

The purpose of this section is to show that Berry's 
cycles studied is this paper are able to entangle N non- 
correlated one-half spin states, while such an effect would 
not be present with linear spin couplings. Since the anal- 
ysis presented here is exploratory, we shall concentrate 
upon the = 4 case which involves already the basic 
features of our entanglement procedure. 

It is convenient to classify the set of the four states 
with respect to the eingenvalues of Sz = J2i^zi- M = 
''^i- The case M = 2 is trivial. The case M — 
1 is the most interesting one for our purpose and will 
be the main subject of this paper. The case M — 
will be discussed briefly at the end of this section. The 
results for negative values of M are readily obtained from 
the positive ones by using the reflection laws introduced 
in section III. The vector space generated by the linear 
combinations of non-correlated spin states with M = 1 
has a dimension four and the natural choice for a basis 
is the set of the four orthogonal states: 



$(1) 

$(2) 
$(3) 
$(4) 



(84) 



We shall assume that the time evolution of the four- 
spin state ^Ji{t) is governed by the time-dependent 
Hamiltonian Hff{t) = iJ(B(t), E(t)), which looks for- 
mally like the quadratic spin Hamiltonian of Eq. ([1} dis- 
cussed extensively in the present paper: 

i?(B(t),E(i)) = U{R{t)) H{B{t),E{t))U-\R{t)), 
H{B{t),E(t)) = -isB{t){Sz/n + X{t){S^/h}2). 

The crucial difference lies in the fact that S is meant to 
he the total spin operator S = X^iLi ^i- 

We shall take as initial state for ^'4(i), one of the 
non-corellated states appearing in equation (j84[) . say, 
^4(0) = ^f^^^. By expanding Hi^{t) in terms of single- 
spin operators one gets a sum of spin-spin interactions: 
cx Si • n Sj • h. where n is given by ri = U{R{t))x. This 
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indicates that Hi^{t) can create spin correlations in the 
state \E'4(t). As a final remark, we stress that, by con- 
struction, H4^{t) is invariant under all permutations of 
the N = A spins. This invariance property will play an 
essential role in our entanglement procedure. 



B. Expansion of a non-correlated four spin state 
into a sum of the eigenstates 



As a first step, we shall construct four angular momen- 
tum eigenstates which are linear combinations of the four 
states Using the rules of adddition of quantum an- 
gular momenta, one finds that the possible eigenvalues of 
, h^S{S + l), correspond to S' = 2 and S' = 1. There is 
a unique way to construct the state S = 2. It is obtained 
by applying the operator S- = Sx — iSy upon the state 
S = 2,M = 2 i.e. *2 2 = I5) «) |i) ® ||) ® One gets 
immediately: 



*2 1 = + + + *^^^). (85) 

The above state is clearly invariant under all permuations 
of the 4-spin states. 

One must now construct three orthogonal states with 
S = M = which will have different symmetry 

under the permutations of the four spins. If one ignores 
the orthogonal condition, it is easy to get three 5 = 1 
states which are linearly independent. The method con- 
sists in taking a linear combination of 2 states in 
such a way as to factor out the two spin singlet state: 
(1^) ®\ — \) — \ — \)® \\))/V2, which is invariant under 
the rotation group. Let us give a typical example of such 
a construction: 



$(1^2) ^ ^^(1) _ ^(2) 



2, 



(86) 



is a good candidate for a S* = M = 1 state. Two 
similar S = M = 1 states $(i.4) are obtained 

from $(^'2) by performing the cyclic permutation (234). 
These three states are linearly independent but a basis 
with orthogonal states would be more convenient. Since 
symmetry will be the basic tool in this subsection, let 
us start with the more symmetric S = M — 1 state : 

1($(1) _ $(2) _^ $(3) _ $(4))_ rpj^g Q^J^gj. g|.g^^gg 

obtained by applying the cyclic permutation (234). One 
arrives then to the following basis for the S — M — 1 
subspace, easily seen to be orthogonal: 



_ $(3) +$(4) _ $(2))^ 
_ $(4) +$(2) _ $(3))^ 



(3) 



The three above states have, different symmetry proper- 
ties under the three permutations (23), (24), (34), which 
are listed below: 



{(23),(24),(34)}vl'\i 
{(23),(24),(34)}M'ii 
{(23),(24),(34)}^Pii 



{-1, 1, -l}*\i, 
{-1,-1, 



The four eigenstates ^ and \1/ 2 1 constitute a complete 
ortho-normal basis for the four spin one-half states with 
M — 1. The expansion of the four non-corrrelated states 
are readily obtained from equations (|84l) and (|85|) 



$(0 

0-3 J 



'j=3 



with 

{1,1,1}, 02, ={-1,-1,1}, 
{1,-1,-1}, a4j = {-1,1,-1}. 



(88) 



(89) 



(87) 



C. The adiabatic evolution of four non-correlated, 
1/2 spin states governed by H{'B{t),E{t)) 

The Hamiltonian H^(t) = i?(B(i), E(t)), as noted be- 
fore, is invariant under either one of the three "pari- 
ties" £(23), £(24), 6(34), associated with the permutations 
(23), (24) and (34). As a consequence, the matrix ele- 
ments: (*\ jv/|i?4(i)l*'\ M') for M = M' = 1 vanish if 
i ^ j . Moreover, the same result holds if ill — ±1, 
M' — ±1 and M 7^ M', since the lowering operator 
Sx — iSy commutes with all the four spin permutations. 
Using the same symmetry arguments, one sees immedi- 
ately that the matrix elements {^2 M\H4{t)\'^\ j^p) van- 
ish whatever the values of i, M, M' . If Ui(t) stands for 
the unitary operator associated with the quantum evolu- 
tion governed by H4{t), the four states U4{t)'^2 m, and 
and U4{t)'^\ j^i will have the same permutation symme- 
tries as their parent states. In conclusion, the four states 
^2M,'I'iA/' behave vis a vis the Hamiltonian H4^[t), as 
if they were associated with isolated spins S and we can 
apply to them the results derived in the previous sections. 

Our Berry's cycle is organized in three steps. At i=0, 
^4(0) = jsB{t)S^. The first step < t <T involves 
the adiabatic ramping of A(t) from to Ao — — 1. In the 
second step T < t < 3T, one proceeds to the rotation 
of E around B by an angle 3tt, while keeping X{t) — Xq. 
For the third step, 3T < t < AT, X{t) makes an adiabatic 
return to its initial value A = 0. During the whole cycle 
the time dependences of X{t) and a{t) are described by 
Blackman functions. The effect of this choice is to tame 
the non-adiabatic oscillating corrections which would be 
generated, if these parameters had discontinuous time- 
derivatives. 

We have performed an exact theoretical analysis of the 
cycle by solving the corresponding Schrodinger equation 
in the rotating frame with jsBT — 25. We have found 
that the adiabatic approximation is working to better 
than 0.1%. It turns out that somewhat accidentally the 
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difference of the dynamical phases for S* = 2, M = 1 and 
S" = 1, M = 1 is close to modulo 27t. Slightly tuning the 
A-rising (-lowering) times can improve the cancellation 
to the 0.1% level. The cyclic evolution of the four spins, 
initially non-correlated, appears then as completely gov- 
erned by the Berry's phases /3(2, 1; Aq) and /3(1, 1; Aq). 



D. The final construction of tiie four-spin 
holonomic entangled state 

At the end of the Berry's cycle, is transformed 
into the following state: 

•fsk^o) = i{exp(z/3(21;Ao))vI/2i + 

exp{^P{ll■Xo))J2a^J^\l}. (90) 



with M=l listed in equation ([84] 



^Bpi^max) = exp(ix) 




(93) 



where x = /3(1 1; Knax)- In order to make more apparent 
the holonomic entanglement resulting from the Berry cy- 
cle, let us rewrite ^gp{Xmax), using the explicit forms of 
the non-correlated states given in Eq. 



B P K'^'f^iax J 



exp(zx) , I 1, 1 1, 1 1, 



1, 



1, 



1, 



,1, 



b>®lo)®b)®|-o>) 



(94) 



Using Eq.® to express the sum Y^^Z^ a,j ^ in terms non-correlated ^spins with M = J2 
of $W and ^ai, one obtains the final expression: 



$W(Ao) = exp(z/3(ll;Ao)) (^<f(i)-i*2i) + 



1 



exp(z/3(21;Ao))«'2i. 



(91) 



(We can verify on this expression that if A/3(Ao) = 
/3(2 1; Ao) - /3(1 1; Aq) = 0, then $^^], coincides with 
up to a phase, as one expects). 

From the energy curves given in Fig.l one sees that the 
region A > should be, in principle avoided, since the 
nearly crossing levels = 2,M = 2 and 5' = 2,M = 1 
could spoil the validity of the adiabatic approximation, 
while the region A < is much more favourable. How- 
ever, in the present context where A and a are the 
only time-dependent parameters, this is only a protec- 
tion against stray magnetic fields orthogonal to the main 
field, since H^^it) has no AM = ±1 matrix elements. The 
two Berry's phases are given explicitly by rather simple 
expressions: 



A similar approach can be used to treat the case of four 

* nii — 0. Among 
the three angular momentum states S = 2, 1,0; Af = 0, 
only the state S — 2,M — has a non- vanishing Berry's 
Phase. As a consequence, one gets an entanglement 
which is weaker than in the above M = 1 case where 
it is maximum. The case of 8 one-half spins (or 8 Qbits) 
with M = 1 is of particular interest since it involves 
four non-vanishing Berry's Phases for S = 1,2,3,4 but it 
will require more powerful Group Theory tools, like the 
Young Operators [30]. However, the basic ideas of the 
entanglement procedure would remain the same. 

The reader might have the impression that the above 
method of entanglement is constrained by the fact that 
Berry's phases depend upon few parameters. However 
one should keep in mind that /3{SM; A) is actually given 
by a Bohm-Aharonov integral involving an Abelian gauge 
field along a closed loop drawn upon the four-dimension 
manifold CP^, which implies a considerable freedom (see 
Eq. (HHj)). 

From a practical point of view entanglement methods 
based on geometric phases are known to present advan- 
tages of robustness: they are resilient to certain types of 
errors in the control of the parameters driving the quan- 



tum evolution [11 1 



/3(21;Ao) = 37r 
/3(ll;Ao) = 37r 




(92) 



In order to obtain the maximum entanglement of 
after a rotation of Sir, it is interesting to choose Aq = 
Xmax = —0.97 which leads to A(3{X„iax) = — tt. Re- 
placing 1 by its expression in terms of <1'^*-' ( Eq. ([85)) ') 
we arrive at the final expression of the four quantum 
entangled states, generated by the Berry's cycles with 
A/3(Ao) = — TT, from any of the four non-correlated states 



We end this section by giving an example where the 
collective spin-spin interaction between a few spins, sim- 
ilar to that described by H]s[{t), has been successively 
implemented. This concerns the case of two or four alkali- 
like ions, without nuclear spin, in an ion-trap. The idea 
is to realize illumination of the N ions (N Qbits) with 
two lasers fields having two different frequencies so that 
the two-photon process, exciting any pair of ions in the 
trap is resonant, but neither of the freq uencies are reso- 
nant with single excitation of an ion [3l[[32j. Realization 
performed for two and four ions in a trap [33| . can be 
generalized to more ions. In the present context it would 
be necessary to match the magnitude of the dc magnetic 
field and the radiation fields. According to the analysis 
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given above, rotating the hnear polarization of the radi- 
ation fields around the B-field should make it possible to 
generate holonomic entanglement of the Qbits. 

VII. SUMMARY AND PERSPECTIVES 
A. Synopsis of the paper 

The purpose of this paper is the theoretical study of 
the Berry's phases, generated in cyclic evolutions of iso- 
lated spins of arbitrary values. The spins are assumed 
to be non-linearly coupled to time-dependent external 
electromagnetic fields (possibly effective) via the super- 
position of a dipole and a quadrupole couplings. Config- 
urations leading to degenerate instantaneous eigenvalues 
are avoided. In other words non-Abelian Berry's phases 
are not considered. We also assume that the two effective 
fields are orthogonal, a mild restriction but with many 
advantages. This implies several discrete symmetries of 
the spin Hamiltonian which simplify considerably the al- 
gebra. For instance, two angular momentum states hav- 
ing different quantum numbers mi and m2 are coupled 
if, and only if, nii — m2 is an even integer. Furthermore, 
for 5 = 1, the geometric space of the Hamiltonian pa- 
rameters is isomorphic to the density matrix space CP^ 
which makes the Berry's and Aharonov-Anandan geo- 
metric phases mathematically identical. 

Using rotation group theory, and the aforementioned 
discrete symmetries, we obtain compact expressions for 
the instantaneous eigenfunctions of the Hamiltonian, la- 
belled by the magnetic quantum number m, valid in the 
limit of small quadrupole coupling. 

We derive an explicit compact expression for the 
Berry's phase /3(to) in terms of the usual Euler angles 
ip, 9, a associated with the trihedron defined by the B, E 
directions, and an extra dimensionless parameter A giv- 
ing the quadrupole to dipole coupling ratio. The result is 
given as a loop integral of an Abelian field along a closed 
circuit drawn upon the parameter space CP^, 

/3(m) — <j) p{m, A)(cos 9 dip + da) — m {dp + da). 

where hp{m,X), the average value of the spin angular 
momentum taken along the B field, is just the gradi- 
ent of the eigenenergy with respect to B. The apparent 
simplicity of the above formula conceals the geometry 
contained in p{m,X). This important feature is clearly 
exhibited in Berry's cycles where A and a are the sole 
varying parameters. The case 5 = 2,to = 0, where the 
Berry's Phase is vanishing for a linear spin Hamiltonian, 
is particularly spectacular: p(0. A) is an odd function of 
A and takes the value 1 when A = 1. This peculiar ge- 
ometry is best illustrated in Fig. [5] where Aa = p{0, A) 
is plotted against the spherical coordinate 9 defined by 
A = -2 cot 9 with < 6* < tt. 

The non adiabatic corrections associated with the Eu- 
ler angles derivatives have been analyzed within the ro- 



tating frame attached to the time-varying fields. The 
Coriolis effect generates an extra magnetic field AB 
which involves a linear combination of the Euler-angle 
time derivatives. The longitudinal component along 
the B field is the only one which survives when a is 
the sole time-dependent Euler angle. The correspond- 
ing Hamiltonian H^^ is devoid of any geometry. As 
a consequence, the phase shift acquired at the end of 
the cycle <?i(B -I- AB) is purely dynamical. The Berry's 
phase is incorporated into the dynamical phase under 
the form of its first-order contribution with respect to 
7] = ~AB///B — (cos 6* ip + a)/{'ysB). The higher-order 
terms give all the non-adiabatic corrections associated 
with a when it is the only varying periodic parameter. 
We have also shown that the subset of these corrections, 
odd under a reversal of 77, cancel exactly for "magic" val- 
ues A = A* (77). This cancellation is implemented in the 
experimental project described in reference [H]. The case 
of the non-adiabatic corrections induced by the trans- 
verse magnetic field ABj_/ B is somewhat more involved 
since it introduces a non-trivial geometry and, as a con- 
sequence, a Berry's phase contribution to be added to 
the one coming from the transverse dynamical phase. An 
explicit evaluation of the complete lowest-order total cor- 
rection A/3_l//3 is found to be oc {'p/jsB)'^ up to a nu- 
merical coefficient computed explicitly for S = 2, m = 0, 
and -1. The results are displayed in Fig. S) 

This discussion hinges upon the implicit assumption 
that the finite value of A is reached by an adiabatic ramp- 
ing governed by H(A(i)) = S^^/h + X{t){S^/h)'^, leading 
to a pure eigen-state V'(A), free of any non-adiabatic pol- 
lution. To analyse this problem, it was convenient to per- 
form the time dependent unitary transformation V{X{t)) 
which transform 'H{X{t)), into a diagonal matrix within 
the angular momentum basis. It follows from the sym- 
metry properties of 7i{X{t)) that V{X{t)) can be taken 
as real and symmetric. The Hamiltonian, which gov- 
erns the spin evolution in the transformed basis, contains 
the non-diagonal term iX{t)V'^ {X{t))-§^V{X{t)). Since 
its time dependence is dominated by the time derivative 
X{t), a linear increase of X{t) would be equivalent to a 
rf pulse with sharp edges leading to the oscillating non- 
adiabatic corrections exhibited in FiglSl The standard 
procedure to smooth them out is to use for X{t) a Black- 
man pulse shape [l^ . The efficiency of the procedure for 
taming the non-adiabatic oscillations is clearly exhibited 
in FiglSJ There is a second assumption implicit in the 
rotating frame analysis, namely that the adiabatic ap- 
proximation be valid for the rotating frame Hamiltonian 
H{t). By solving numerically the Shrodinger equation 
describing the whole cycle, we have verified (see ref. [1] 
for explicit numerical results) that, indeed, the adiabatic 
approximation works beautifully provided one also uses 
a Blackman-pulse shape for the angular speed, a{t). 

In the last section, we propose a procedure to en- 
tangle a system of N non-correlated one-half spins (or 
N Qbits). It involves Berry's cycles generated by an 
Hamiltonian formally identical to one given in equation 
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([T|), but with an important change: S stands now for 
the total spin operator S = J2i=i^i- We have used a 
method which is well adapted, to the simple configura- 
tion of four one-half spins. In view of possible extensions 
to N spins with A'^ > 4, we would like to rephrase our 
approach within a more general framework. A configu- 
ration of N non-correlated ^-spins can be described by 
a factorizable spin tensor of order N. It constitutes a 
reducible representation of the S{U2) group. There is 
a general procedure to decompose this tensor into irre- 
ducible tensors, associated with eigenstates of and 
Sz- With the exception of the trivial case N — 2, the 
angular momentum state vPsm can appear several times 
in the decomposition. According to a known theorem 
of Group Theory [s^ , all the states 4'^ j^j have different 
symmetry properties under the permutation of the N 
spins and they can be organized into an orthogonal ba- 
sis, (^smI^S'M') ~ Ss^s'^ALM'Sij- Let us now consider 
the set of Hamiltonians H(S) which are given by a non- 
linear series expansion with respect to the spin compo- 
nents Sx.y.z- By construction, 'H(S) is invariant upon any 
permutation of the TV spins. As a consequence, 'H(S) is 
diagonal within the basis ^"5^: {^^g ^\'H{S,t)\^-'g, j^.^,) = 
i5s,S'i5ij(*SA/l^(^'^)l*s A/')- The above equation ex- 
presses the simple fact that 'H(S) is acting upon the states 
M if they were isolated spins S. Any initial state of 
TV non-correlated spins can then be written as |vE'Ar(i — 
0)) = s'^s'mI^sm)- One performs now an adiabatic 
Berry's cycle along a closed loop in the parameter space 
such that ?^(S,0) =n{S,T) =7sS-B. At the end of the 
cycle the N spin system is given by the following state 
vector: |*Ar(r)) = Ei. s «SM exp M)|M'^_ ^). 

With a suitable choice of the Berry's cycle, we have 
shown in the particular case N — 4, M — 1 that the final 
state is endowed with a maximal entanglement. Thus 
extension to higher values of N is worth pursuing, re- 
membering that according to quantum computing ex- 
perts: "entanglement, as with most good things, it is 
best consumed in moderation" [34i] . 



splitting, n the Rabi frequency relative to the transition 
55'i/2 — > 5P1/2 and e the light polarization. However 
there is a certain price to be paid arising from the in- 
stability of the "dressed" ground state depicted by its 
decay rate Tdec = -^;r^^5Pi/2 (where T^p^^^ is the spon- 
taneous decay rate of the excited state) . Although, with 
our choice of detunings, this effect does not affect the 
Berry phase value, it can perturb seriously the measure- 
ment process. One is clearly facing a difficult optimiza- 
tion problem, if one wants also to keep the non-adiabatic 
correction below the 0.1% level. We give here some fea- 
tures of the solution described in ref. With an ex- 
ternal magnetic field of 1 mG, one can adjust {i.e. 
the light beam intensity) in order to get a quadrupole to 
dipole couplings ratio close to 1. The off-diagonal den- 
sity matrix element holding the Berry Phase, decays to 
half its initial value. But adjusting the interferometer 
parameters can compensate for this effect and keep the 
interference contrast close to its maximum. 

On the above example the role of both the hyperfine in- 
teraction, and the instability of the atomic excited states 
is clearly exhibited. This gives a clear illustration of 
the atomic internal structure contribution to the spin 
dynamics. Although these effects upon Berry's cycles 
can be accounted for by choosing appropriate values of 
the effective B, E fields, they lead to severe experimental 
constraints. It looks possible to satisfy these constraints 
with *^Rb atoms but with light alkalis this is far from 
obvious Q. On the other hand, the external degrees of 
freedom, i.e. the atomic nuclei coordinates have been 
supposed fixed physical quantities, an assumption which 
would obviously fail for ultra cold atoms belonging to a 
fermionic or bosonic quantum gas. In such a case, these 
coordinates become truly quantum variables which ap- 
pear explicitly in the density matrix of the system. The 
problem of Berry's phase generated by an adiabatic cycle 
of the coupled external fields is well beyond the scope of 
the present paper. 
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B. Experimental perspectives 

In a separate work, we have explored the possibility 
of cold alkali atoms in their ground state to measure 
Berry's phase with atomic interferometers. The spin 
operator is then identified with the total angular mo- 
mentum F = s -I- I. The ac-Stark shift induced by a 
linearly polarized light beam tuned off-resonance of one 
resonance line, can induce the quadratic interaction if 
one accepts a few experimental compromises. The can- 
didate for our spin system is the ground state hyperfine 
(hf) level 55i/2, F = 2 of s'^Rb. We have found judi- 
cious to tune the laser frequency midway between the 
two lines SSi/a, F = 2 ^ 5Pi/2, F = 1, and F = 2 . 
The effective quadrupole coupling takes then the simple 
form: Hq = 



(F • e)^ where AWp is the 5Pi/2 hf 



Appendix. Some useful formulas and algebraic 
results for S=3 and S=4. 
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det(^odd(3, A) - a; 1) = a;^ - Ux^X + 49a;^A^ - lOa;^ - 36a;A^ + 102a;A - 216A2 + 9, 



(97) 



dct(?^e,;en(3, A) - a; 1) = a;^ - 14a;2A + AQxX"^ - 4x - 36X^ + 24A. 



(98) 



/ 2A + 4 V7X \ 

VlX 8A + 2 S^lA 

S^lA lOA S^fA 

S^fA 8A-2 VlX 

\ V7X 2X-A J 



(99) 



nodd{A,x) = 



/llA+3 3^ \ 

+ 1 5A 
5A i|A _ 1 iH^ 
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(100) 



det(We«en(4,A) -a; 1) = a;^ - 30a;^A + 273a;^A2 - 20a;^ - 820a;^A^ + 472a;^A, (101) 
det(^odd(4, A) - X 1) = x'' - 30a;3A + 273a;2A2 - lOa;^ - 820a:A^ + 182a;A + 576A'^ - 712A2 + 9. (102) 



[1] M.V. Berry, Proc. R. Lond. A392, 45 (1984). 
[2] B. Simon, Phys. Rev. Lett. 51, 2167 (1983). [18 
[3] Y. Aharonov and J. Anandan, Phys. Rev. Lett. 58, 1593 [19 
(1987). 

[4] N. F. Ramsey, Phys. Rev. 76, 996L. and Molecular Beams [20 
(Clarendon, Oxford, UK, 1956), pp. 124-133 (1949). [21 
[5] M.A. Bouchiat and C. Bouchiat, to be pubUshcd. [22 
[6] A. Shapcrc and F. Wilczek, Geometric Berry's Phases in [23 
Physics, Advanced Series in Mathematical Physics, Vol. 
5 (World Scientific, Singapore) (1989). [24 
J. Anandan, J. Christian and K. Wanelik, Resource letter 
GPP-1 Geometric phases in physics, Am. J. Phys. 65, [25 
180 (1997). [26 
M. Berry, Sci. Am. 259, 46 (1988). 
[9] Holstein, J. Am. Phys. 57, 1079 (1989). [27 
[10] J. Zwanziger, M. Koenig and A. Pines, Ann. Rev. Phys. [28 
Chcm. 41, 601 (1990). [29 
E. Sjoqvist, Physicsl, 35 (2008). [30 
J. A. Jones, V. Vcdral, A. Ekert, and G. Castagnoli, Na- [31 
turc 403, 869 (2000). 
[13] E.D. Commins, Am. J. Phys. 59, 1077 (1991). [32 
[14] M. Pendlebury et al, Phys. Rev. A 70, 032102 (2004). 
[15] M.R. Tarbutt, J.J. Hudson, B.E. Sauer, E.A. Hinds, [33 
Faraday Discuss. 142, 37 (2009). [34 
[16] T. Bergmann, T. Gasenzer, O. Nachtmann, Eur. Phys. 

J. D 45, 197 and 211 (2007). 
[17] DeKieviet et al., 2010, e-print arXiv:physics.atom- 



[7] 



[8] 



[11] 
[12] 



ph/1003.0622. 

R. Tycko, Phys. Rev. Lett. 58, 2281 (1987). 

J.W. Zwanziger, M. Koenig and A. Pines, Phys. Rev. A 

42, 3107 (1990). 

F. Wilczek and A. Zee, Phys. Rev. Lett., 58, 2111 (1984). 

E. Kiritsis, Commun. Math. Phys. Ill, 417 (1987). 

A. Zee, Phys. Rev. A 38, 1 (1988). 

C. Cohen Tannoudji and J. Dupont-Roc, Phys. Rev. A 
5, 968 (1972). 

C. Bouchiat and G. W. Gibbons, J. Phys. France 49, 187 
(1988). 

C. Bouchiat, J. Phys. France 50, 1041 (1989). 

L.I Schifl', Quantum Mechanics, 2nd erf., Sec. 31 

(MacGraw-HiU, New- York, 1955). 

F. J. Harris, Proc. IEEE, 66, 51 (1978). 

M.V. Berry, Proc. R. Soc. Lond. A 414, 31 (1987). 

Y. Aharonov and D. Bohm, Phys. Rev. 115, 485 (1959). 

M. Hamermesh, 1962, Group Theory, (Pcrgamon). 

K. Molmer and A. Sorensen, Phys. Rev. Lett. 82, 1835 

(1999). 

A. Sorensen and A. Molmer, Phys. Rev. Lett. 82, 1971 
(1999). 

C. A. Sackett et al.. Nature 404, 256 (2000). 

D. Gross, S.T. Flammia and J. Eisert, Phys. Rev. Lett. 
102, 190501 (2009). 



